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Abstract 


Pulse-like acoustic signals are transmitted from an acoustic source near Oahu to seven 
receivers off the west coast of the United States for a 124-day period in 1988. Acoustic 
travel-time oscillations are observed in the received signal at periods between 15 and 23 
hours, which are caused by barotropic (or first or second mode baroclinic) fluctuations 
in the ocean. It is shown that these fluctuations cannot be local processes isolated to 
either the source or to the receivers. It is further shown that resonant barotropic gravity 
wave modes (Platzman et al., 1981) are not consistent with the data. The cause of these 
fluctuations remains unresolved, but the data and other oceanographic measurements put 
many constraints on the process causing these fluctuations. 
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Chapter 1 


Introduction 


1.1 Motivation 


Acoustic tomography has been successfully used to analyze mesoscale (order 100-1000km) 
features of the ocean (Spiesberger et al., 1983; DeFerrari and Nguyen, 1986; and others). 
For these experiments, other synoptic methods of extracting information about mesoscale 
features are available (i.e. XBT’s, CTD's, and current meters). Yet when considering 
synoptic measurements on a basin-scale, the cost of deploying such equipment becomes 
prohibitive. In the past, analytic and numerical methods have been relied upon to provide 
dynamic models of these features. An acoustic tomography array can, at reasonable 
cost, be used to obtain synoptic measurements of barotropic, basin-scale oceanographic 
phenomena (Spiesberger et al., 1989b). 

A four month basin-scale tomography experiment was performed in 1988 in the north- 
east Pacific where a bottom-mounted source near Oahu transmitted a continuous (13312) 


acoustic signal to seven bottom-mounted receivers off the west coast of the United States. 





Signal processing converted the signal at the receiver into an equivalent signal which would 
be received if the source transmitted a pulse-like signal of 60 ms duration at 184.2 s inter- 
vals. Changes in acoustic travel-time are nieasured by estimating the phase change of the 
received signal of adjacent records. This experiment follows up two earlier experiments 
(in 1983 and 1987) in which a source near Oahu transmitted to one receiver off the coast 
of northern California for 5 and 21 days, respectively. 

These earlier experiments yield the following results. Travel-time oscillations at tidal 
frequencies are observed, and are caused by barotropic tidal currents. In addition, sig- 
nificant oscillations are observed at may non-tidal frequencies. Specifically, at periods 
between 15 and 23 hours, oscillations are observed having rms variations of about 10 ms 
which vary in amplitude and frequency on about a 10-day time scale, indicating that they 
may be intermittently excited. Since data were only available at one receiver, it has not 
been possible to match these oscillations to an oceanographic process (Spiesberger et al., 
1989b) 

With results from the 1988 experiment, presented in this thesis, we examine four 
hypotheses in a search for the oceanographic source of these oscillations. The hypotheses 


are: 


The ocean fluctuates only near one receiver. 


The ocean fluctuates independently at each receiver. 


The ocean fluctuates only near the source. 


The ocean fluctuates at basin-scales as a result of resonant barotropic gravity wave 


modes. 





We show each of these hypotheses to be unlikely, and are not able to provide a theory 
for the cause of the observed travel-time oscillations which is consistent with both the 


tomographic data and other measurements of the ocean. 


1.2 Overivew 


The thesis is organized as follows. Chapter 2 contains a brief discussion of acoustic 
propagation in the ocean as it relates to acoustic tomography. Chapter 3 discusses the 
experiment and signal processing. Chapter 4 examines each of the four hypotheses pro- 
posed to explain the cause of the acoustic travel-time oscillations. Chapter 5 concludes 


with a discussion of the constraints placed on the problem by both the tomographic data 


and other measurements, and suggests other avenues of research for this problem. 


10 





Chapter 2 


Acoustic Propagation 


l in water, which is about four and one-half times that 


Sound travels at about 1500 m-s” 
of air. Sound speed increases with increasing temperature, increasing pressure, or, to a 
much smaller extent, increasing salinity. At middle and equatorial latitudes, the speed of 
sound initially decreases with depth because of cooling, and then increases in the deep 
ocean due to great pressure, resulting in a sound speed minimum at about 1 km depth 
(Figure 2-1). Application of Snell's law of refraction shows that sound is refracted towards 
regions of lowest sound speed (Officer, 1958). The ocean thus has a sound channel (Ewing 
and Worzel, 1948) which focuses acoustic energy away from the surface and bottom where 
significant attenuation and scattering of the signal would otherwise occur (Figure 2-1). 
This fact, coupled with remarkable water transparency to sound at frequencies below 
about 300 Hz (Urick, 1983), allows detection of sound propagated over thousands of 
kilometers. 


At frequencies above about 100 Hz, propagation of sound can be modeled with ray 


theory, an approximate solution of the acoustic wave equation, since the acoustic wave- 
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Figure 2-1: Vertical profile of sound speed in the northeast Pacific half-way between Oahu 
and California (left). Twelve acoustic ray paths originating from a source on the sound 
channel axis for a 200 km section of ocean (Tight). 


length (about 15 m at 100 Hz) is much less than the vertical scale of the channel (about 1 
km). Ray theory predicts that an acoustic pulse emitted from a submerged source travels 
to a distant receiver along ray paths (multipaths) which have different travel-times in ac- 
cordance with Snels law of refraction (Officer, 1958). Each ray path samples a different 
vertical slice of the ocean, based on the departure angle (from the horizontal) of the ray 
from the source. Thus the travel-times of different multipaths can be used in an inverse 
problem to solve for baroclinic fluctuations in the ocean. Alternatively, travel-times from 
all multipaths can be combined in an inverse problem to solve for barotropic fluctuations 


in the ocean. 


Yo first order, fluctuations in acoustic travel-time can be computed by assuming that 


the ray path does not move in the presence of fluctuations (Munk and Wunsch, 1979). 


th 


Perturbations of acoustic travel-time, 67;, for the 2°" ray path, [;, are primarily related 


EZ 





to in-situ fluctuations of temperature, 60, hydrostatic pressure, pgrj, and current, à, ac- 
cording to, 


QUT) c n( s, t) 


ü(s,t) - ds 


6 {t — — E LS - / = , 
Tilt) af. e TE r, EED 


(2.1) 


where the i” ray path has coordinates parameterized by 5 (with differential elements ds 
and ds), and tis time. The reference sound speed is c,(5,t), and that portion of the 
sea-surface displacement related to hydrostatic fluctuations 1s 9. The water density is p 


and the gravitational acceleration is g. The constants are, 
a 42 110 Hinc spe 


A a 


(Spiesberger et al., 1989b). 

Suppose a ray path travels a length L (km) in the presence of a uniform oceanic fluc- 
tuation where current (parallel to the ray path) is u (m/s), or where the temperature 
fluctuation is 60 (°C), or where hydrostatic pressure changes are caused by sea-surface 
displacement 7 (m). Then if currents, temperature, or sea-surface oscillations are indi- 


vidually responsible, the travel-time oscillation (s) is approximately, 


ór = —444 x 10 ^ Lu (2.2) 
os Loe (2:3) 
ór& —7.55 x 10-9; (2.4) 


(Spiesberger et al., 1989b). 
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Chapter 3 


The Experiment and Signal 


Processing 


3.1 Experimental Layout 


A continuous, phase-modulated signal (center frequency 133 Hz, bandwidth 17 Hz) is 
transmitted between 6 January and 9 May 1988 (124 days) from a bottom-mounted source 
off the north coast of Oahu (183 m depth) to seven bottom-mounted receivers off the west 
coast of the United States (Figure 3-1). The receivers are about 3000 to 4000 km from 
the source and have acoustic travel-times of about 30 to 40 min respectively. The source 
is within 30 m of the sources used for the 1983 and 1987 experiments (Spiesberger et al., 


1989b). 
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Figure 3-1: The 1988 acoustic tomography experiment consists of a bottom-mounted 
source (S) located near Kaneohe Bay, Oahu and seven bottom-mounted receivers (R1 to 
R7) whose approximate positions are shown. The distances between the source and the 
receivers are about 3000 to 4000 km. 





3.2 Signal Processing 


In order to raise signal-to-noise ratio above a detectable threshold, the transinitted signal 
is phase-modulated every eight cycles using a 511-digit pseudo-random code. A matched 
filter at each receiver compresses the energy of the 30.7 s cycle into that of an equivalent 
transmitted pulse of 60 ms duration. Six such cycles are averaged together, yielding travel- 
time records at 184.2 s intervals, further increasing signal-to-noise ratio. This processing 
is the same as that used in prior experiments (Spiesberger et al., 1989a; Metzger, 1983). 
At each receiver, we measure changes in acoustic phase, 6@, at intervals of 184.2 s 
using a maximuim-likelihood estimator which computes the phase of the complex cross- 
correlation of adjacent travel-time records (Spiesberger et al., 1989a). This estimator 
is optimally designed to detect barotropic fluctuations as well as first or second mode 
baroclinic fluctuations, and is designed to attenuate signals from higher vertical modes. 
If the acoustic travel-time fluctuates at period T', then the change in travel-time is 


related to the change in acoustic phase by, 


Too 
=. m 3.1 
d (UNT on 


where the source frequency, w., equals 27-133 rad-s7! and the time between travel-time 
records, At, equals 184.2 s (Spiesberger et al., 1989a). Travel-time changes are related to 
oceanographic fluctuations according to Equation 2.1. 

We are interested in oceanographic fluctuations other than internal waves of higher 
mode number (> 2). The acoustic signals associated with these internal waves are con- 
sidered to be a portion of the “noise.” The precision of the travel-time estimates (135ts 
at At — 184.2 s) is limited by the higher mode internal waves rather than by clock errors 
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at the source and receiver or by acoustic noise (Spiesberger et al., 1989a). 
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Chapter 4 


Hypotheses 


We consider four hypotheses to explain measured fluctuation m acoustic travel-time at 


periods between the semi-diurnal and diurnal tides. They are: 


1. The ocean fluctuates only near one receiver (receiver R3 of Figure 3-1). Previous 
papers only analyzed data from R3 (Spiesberger et al., 1989a,b). In section 4.1, we 


find this hypothesis to be false. 


2. The ocean undergoes independent fluctuations near each receiver. In section 4.2, we 


discuss evidence which weighs against this hypothesis. 


3. The ocean fluctuates only near the acoustic source. In section 4.3, we discuss evi- 


dence which disagrees with this hypothesis. 


4. The ocean fluctuates at basin-scales. In section 4.4, we present evidence that the 


oscillations are not resonant barotropic gravity wave modes of the world’s oceans. 


— 
Qo 





4.1 Hypothesis 1: Ocean Fluctuations Localized to R3 


All seven receivers exhibit proininent oscillations in acoustic travel-time between the semi- 
diurnal and diurnal tidal periods (Figures 4-1 and 4-2). Therefore this hypothesis is not 


true. 


4.2 Hypothesis 2: Independent Ocean Fluctuations Near 


Each Receiver 


We first check if acoustic travel-time oscillations of the same frequency are simultaneously 
observed at all receivers. The 124-day experiment is divided into twelve consecutive 
10-day segments. Periodograms are computed for each receiver which collected data 
during the segment. We find that the oscillations between periods of 15 and 23 hours are 
simultaneously observed by most receivers (Figure 4-3). For example, between yeardays 
70 and 79, all seven receivers exhibit an oscillation at periods between 16.1 and 17.5 
hours and exhibit an oscillation at periods between 19.0 and 20.9 hours (Figure 4-3). We 
conclude that most receivers simultaneously observe travel-time oscillations at the same 


frequency. 


We believe the above data strongly weighs against hypothesis 2. Further evidence 
against this hypothesis is presented next. 
The rms variations of acoustic travel-time are about 10 ms for each peak observed 


in the 10-day periodograms at R1 to R7 (Figure 4-3). In 1983 and 1987, acoustic os- 
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Figure 4-1: Power spectral densities of a 42-day time series (yeardays 41 to 83, 1988) 
of phase differences at each receiver. A Bartlett window is applied to the data. The 
mean spectral level due to random walk noise (Spiesberger et al., 1989a) is at or below 
1 x 107°cycles?/cpd. The vertical solid lines mark prominent tidal periods. The vertical 
dashed lines mark periods of all resonant barotropic gravity modes with periods between 
15 and 23 hours (Platzman et al., 1981). 
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Figure 4-2: RMS travel-times (ms) are shown for the peaks in Figure 4-1 with periods 
between 15 and 23 hours (vertical axis). The arrows (right axis) mark the periods of all 
resonant barotropic gravity modes with periods between 15 and 23 hours (Platzman et 


al., 1981). 
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Figure 4-3: The number of receivers exhibiting spectral peaks within a given range of 
periods is shown for each 10-day segment of the experiment. The number of receivers for 
which data are available for each segment is noted in parenthesis in the left column. In 
each case, the spectra are computed using a Bartlett data window. The periods listed 
(top) are the eigenfrequencies of the periodogram. The spacing between eigenfrequencies 
is approximately half the resolution bandwidth of the periodogram (Harris, 1978). 
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cillations near a 15 hour period had similar size (Spiesberger et al., 1989a,b). Applying 
Equations 2.2 and 2.3, we calculate how large a barotropic current or how large a temper- 
ature perturbation must be to yield an rms variation of 10 ins. If the ocean fluctuations 
are confined within a 50 km region of each receiver, a uniform current of 0.6 m-s7! or a 
uniform temperature perturbation of 0.132? C is required. Magnitudes of the perturbation 
in current and temperature are inversely proportional to the horizontal scale of the fluctu- 
ation. If the fluctuation is confined to 25 km, then the associated current or temperature 
perturbation would be twice the above values. The magnitudes of these perturbations are 
implausible. 

We therefore find it unlikely that independent ocean fluctuations near each receiver 


give rise to the observed oscillations in acoustic travel-time at periods between the semi- 


diurnal and diurnal tides. 


4.3 Hypothesis 3: Oceanic Fluctuation Localized Near the 


Acoustic Source 


We present two results which indicate that the acoustic oscillations with periods between 
15 and 23 hours are not caused by ocean fluctuations confined to the region of the acoustic 


source. These results are presented separately in subsections 4.3.1 and 4.3.2 below. 


4.3.1 Island Trapped Waves 


A model for the trapping of super-inertial internal waves around islands was proposed 


by Wunsch (1972) to explain thermal fluctuations near Bermuda. In this model the 
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ocean (constant depth, rigid lid), containing a cylindrical island, is realistically stratified 
allowing for free internal waves. Internal wave scattering exhibits pseudo-resonauces if tlie 
period of the incident wave matches the free wave travel-time around the cylinder (island ). 
When this model is applied to Oahu (Appendix A), no significant pseudo-resonances are 
observed for wave periods between 10 and 25 hours. 

A second mechanism for the generation of island trapped waves, refraction (resulting 
in “trapped-leaky” waves), is not likely to occur around Oahu due to the narrow shelf 
width and complicated geometry of the island (Luther, 1985; Longuet-Higgins, 1967; 
Summerfield, 1972). Thus, available theoretical models fail to provide a mechanism for 
the generation of island trapped waves around Oahu. 

Even though we find no theoretical justification for island trapped waves, these theories 
are only approximations which contain some questionable assumptions, and such waves 
may exist. For example, Luther (1985) has analyzed 10- and 14-year records of sea-level 
at two stations on opposite sides of Oahu (Mokuoloe and Honolulu). His analysis reveals 
peaks in the sea level spectra from each station at periods of 17.0 and 20.0 hours with 


amplitudes of about 0.5 cm. Cross-spectral phase relationships suggest that these spectral 


1 


peaks could result from island trapped waves with horizontal wavelengths equal to 1 or 5 


the circumference of Oahu. Patzert and Wyrtki (1974) have observed clockwise currents of 
about 5 cm/s around Oahu, providing further support for the existence of island trapped 
waves. 

If there is an island trapped wave of the type hypothesized by Luther, this is not 
the wave observed in the tomographic measurements. Consider the transmission paths 


to R3 and R4 which leave the source at bearing angles differing by about 10° (Figure 
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4-4). If a 15.5 hour period island trapped wave has a horizontal wavelength of 5 the 
island circumference, then the predicted phase lag between acoustic signals at R3 and 
R4 is 5° (Figure 4-4). Cross-spectral phases at 15.5 hour period (at R3 and R4) have 
values of 157? + 15° (Figure 4-5 and Appendix B). Calculations of cross-spectral phases 
at other receiver combinations and periods yield results which are also inconsistent with 
the island trapped wave hypothesis of low horizontal mode number (Figure 4-5). If the 
horizontal mode number of the island trapped wave is 64, then the predicted phase lag at 
R3 and R4 would match the observed phase lag. We do not believe that realistic coastlines 
and bathymetric features could support such an island trapped wave. The tomographic 


measurements are therefore inconsistent with the island trapped wave hypothesis. 


4.3.2 Other Near-Source Fluctuations 


Unrealistically large fluctuations of current and temperature are required to yield acoustic 
travel-time oscillations of 10 ms rms (Section 4.2). We conclude that the tomographic 


measurements are inconsistent with hypothesis 3. 


4.4 Hypothesis 4: Resonant Barotropic Gravity Wave Modes 


We investigate the hypothesis that the acoustic oscillations are caused by resonant barotropic 


gravity wave modes. 





tom 


to R4 





N WAVE PATH 


Figure 4-4: The relationship is shown between the phase lag, n®, of an island trapped 
wave (horizontal mode n) at the wave's intersection with two ray paths and the difference 
in bearing angle of the ray paths, O. The tomographic signal generated by an island 
trapped wave will have a phase lag between the two receivers equal to about no, where 
nd — no. The distances r; and r2 are about 40 and 10 km respectively, and O is about 
10°. The angles in this figure are enlarged for clarity. 
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Figure 4-5: Phase lags (degrees) are shown for periods at which coherency peaks occur 
in the cross-spectral estimates of 42-day time series of phase differences at each receiver 
referenced to R3. Each time series has the same reference time and runs from yearday 41 


to 83. 1988. The error bars show the range of periods over which the phase shift remains 
: 


within 15º of its value at the coherency peak (Appendix B). The arrows (right axis) 
mark the periods of all resonant barotropic gravity modes with periods between 15 and 


23 hours (Platzman et al., 1981). 
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4.4.1 Resonant Shallow Water Gravity Modes 


Platzman (1978) and Platzman et al. (1981) developed a model identifying normal modes 
of the world ocean. In this model, the primitive equations are applied to a finite element 
discretization of the world ocean (about 900 km between nodes) yielding eigenfrequencies 
of resonant barotropic vorticity and gravity modes. Corresponding to each eigenfrequency 
is a set of eigenvalues describing the sea level and velocity structure of each mode. Platz- 
man et al. (1981) find eight gravity modes at periods between 15 and 23 hours (Figures 
4-1, 4-2, and 4-5), two of which (15.5 and 21.2 hours) are dominant in the Pacific Ocean. 

In our study, we use the same programs developed by Platzman to identify resonant, 
barotropic gravity modes for the Pacific Ocean only (Miller et al., 1988). Using the 
same grid as Platzman (with minor variations at boundary nodes), we find three resonant 
modes at periods between 15 and 23 hours (15.06, 18.46, and 21.66 hours). Considering 
the altered boundary nodes, this 1s consistent with the Platzman’s results. 

For each mode, we compute the predicted acoustic travel-time oscillation, at receiver 


?, using Equation 2.1 and obtain results of the form, 


ór;(t) 2 aajsin(-Eat 4 ój - C), (4.1) 


where a; and $; are the amplitude and phase of the acoustic travel-time computed from 
the model, o is frequency and t is time. Since the model magnitude and phase of the sea- 
surface and velocity are based on arbitrary references, scaling factors a and C are required 
when comparing the predicted and observed acoustic travel-time oscillations. For each 
mode, ¢ is chosen such that the phase of the predicted acoustic travel-time at R3 (the 


reference receiver) is 0°, and a is chosen such that for the amplitude of the predicted 
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acoustic travel-time matches that of R3. We estimate the accuracy of our model travel 
times to be about 20% in magnitude and +15° in phase. 

A comparison of the observed and predicted acoustic travel-times at each receiver are 
shown in Tables 4.1, 4.2, and 4.3 for each of the three modes. The model phase at each 
mode deviates significantly from the observed phase. For example, at the 15.06 hour 
mode there is a consistent —160? discrepancy between the model and the observations. 
We conclude that the observed travel-time oscillations are not caused by a single mode. 

We investigate the case where two modes occur simultaneously at each modal fre- 
quency (clockwise and counterclockwise). We find no combination of these modes at any 
of the three modal frequencies yield the observed phase and magnitude of acoustic travel- 
time oscillations at each receiver (Appendix C). We therefore conclude that the two-mode 


hypothesis is not consistent with the observations. 


4.4.2 Comparison With Other Measurements 


For each resonant barotropic shallow water gravity mode modeled in section 4.4.1, we 
observe that in order to yield an averaged (42-day) rms acoustic travel-time oscillation of 
about 3 ms, sea surface oscillations of about 5 cm are required. However, other obser- 
vations of the sea-surface at these periods and record lengths exhibit magnitudes which 
are about ten times less (Table 4.4). Errors from boundary condition approximation and 
discretization cannot account for this difference. Based on this and on the failure of the 
model to fit the data, we conclude that the tomographic oscillations observed are not the 


result of resonant barotropic gravity wave modes of the Pacific Ocean. 
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15.06 Hour Mode 
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15.06 Hour Mode 


Positive Frequency Negative Frequency 


Receiver | Observed Phase | Model Phase Deviation | Model Phase Deviation 


Lag (deg.) Lag (deg.) (deg.) Lag (deg.) (deg.) 


-161 | -157 
-169 -162 


0 0 


70 -140 
91/116 | 177/152 
70/-86 2167/21] 

19167 68/-144 





Table 4.1: Comparison of acoustic travel-time amplitudes (top) and phases (referenced to 
R3) (bottom) predicted by the resonant barotropic gravity wave model and tlıose observed 
in the tomographic data at the 15.06 hour mode (Miller et al., 1988). Phase deviations are 
determined by subtracting the observed phase lag from the model phase lag. Observed 
phases are modulated by 4360? if necessary to obtain deviations between X180? (a —175? 
deviation is only 10? different from a 4-175? deviation). Positive and negative frequency 
indicates clockwise and counterclockwise wave propagation. Phase lags are accurate to 
+15° (Appendix B). When two auto- or cross-spectral peaks occur near the predicted 
mode period, both are tabulated. 
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18.46 Hour Mode 


Observed Model 
Receiver | Amplitude (ms) | Amplitude (ms) 


18.46 Hour Mode 
Positive Frequency Negative Frequency 


Receiver | Observed Phase | Model Phase Deviation | Model Phase Deviation 


Lag (deg.) Lag (deg.) (deg.) Lag (deg.) (deg.) 


03 131 115 
150/115 -141/-106 | -159/-124 
0 0 0 
2120.82 111/-141 129193 
-6/-154 -9/139 21/169 
126/8 -144/-26 -108/10 
69 -65 T 





Table 4.2: Same as Table 4.1, except comparisons are made at the predicted 18.46 hour 
mode. 
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Observed 
Amplitude (ms) 


Receiver 





21.66 Hour Mode 


Receiver | Observed Phase 


Lag (deg.) 


Model Phase 
Lag (deg.) (deg.) 
63/174 67/178 
be al -60/-73 
0 0 
35/180 -30/-175 


93 -39 
-93 63 


77/-107 -85/99 


Model 
Amplitude (ms) 


Deviation 








Positive Frequency Negative Frequency 
Model Phase 


Deviation 
Lag (deg.) (deg.) 
-59/-170 
-56/-69 
0 
-40/175 
-67 
43 
-69/115 


Table 4.3: Same as Table 4.1, except comparisons are made at the predicted 21.66 hour 


mode. 
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Figure 4-6: Power spectral densities of 42-day time series of sea level at La Jolla, California 
(top) and Mokuoloe, Oahu (bottom). The La Jolla sea level data is from yeardays 41 to 
83, 1981 and the Mokuoloe sea level data is from yeardays 19 to 61, 1988. A Bartlett 
window is applied to the data. The labeled peaks at periods between 15 and 23 hours are 
not sidelobes of the tides. 
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Source (Reference) 


Oahu sea-level 10-year record (Luther, 1985) 


Oahu sea-level 12-day record (Figure 4-6) 


La Jolla, Calif. sea-level 42-day record 
(Figure 4-6) 

Bottom pressure gage (300 km east of Oahu) 
60-day record (Filloux, 1981) 

Bottom pressure gage (1150 km west of Calif.) 
7-day record (Filloux, 1971) 





Inter-tidal RMS Amplitude 
Peaks Observed? Oscillation 


0.5 cm 


0.5 cm 


0.2 cm 


0.1 cm 


Table 4.4: A summary showing amplitudes of sea level oscillations above background 


spectral levels at periods betwcen 15 and 23 hours. 
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Chapter 5 


Discussion 


5.1 Summary of Results 


We have shown in the previous chapter how the tomographic data provides several con- 
straints on the problem of identifying a cause for the observed travel-time oscillations. 
Other experiments in the Pacific Ocean add further constraints. At this point, having 
eliminated the more obvious processes which could cause these oscillations, it is useful to 
summarize these constraints and perhaps obtain a new perspective on the nature of an 


oceanographic process required to generate the observed results. Such a process must: 
e have about 0.5 cm oscillations in sea level. 


e generate, through current or temperature perturbations, acoustic travel-times of 
about 10 ms rms over 10-day time periods and about 3 ms rms over 42-day time 


periods. 





e result in large phase lags in acoustic travel-time between different receivers consistent 


with tlıe data. 


e be intermittent (averaging over time reduces the magnitude of the oscillations). 


e occur simultaneously at periods between 15 and 23 hours at similar magnitudes. 


The requirement for small sea level variations brings us back to the concept of internal 
waves. Because of the signal processing, our data is insensitive to higher order internal 
wave modes, and we have previously ignored such processes. ‘There is a possibility that 
internal waves at the turning latitude (where wave frequency equals the inertial frequency) 
may intensify and give rise to temperature perturbations sufficient to yield the observed 
travel-time magnitudes (M. Briscoe, personal communication ). 

The inertial period in the region of the experiment ranges between 32.7 hours (at 
Oahu) and about 15.7 hours (at R1). In order for travel-time fluctuations to occur at a 
receiver due to this process, the ray path must cross the turning latitude for a particular 
period. For example, at R5 and R6 (inertial period about 24 hours), we expect to observe 
fluctuations due to this process only between 24 and 32.7 hour periods. Yet we observe 
prominent oscillations at these receivers at 15 to 23 hour periods. Any test of this theory 
must account for these oscillations. 

For completeness, we note in the literature that several other theories explaining the 


observed travel-time oscillations have been examined and rejected (Spiesberger et al., 


1989b). 


5.2 Conclusion 


We must confess that, while the results presented here are significant, they are unsatisfying 
in that we are unable to find an oceanographic process which explains the acoustic travel- 
time oscillations. Had it not been for other work demonstrating agreement between the 
tomographic data and models of barotropic tides (Headrick, 1990), we may have even 
questioned the validity of the data itself. In one sense, this study highlights the concept 
that while it is easy to design an experiment to prove or disprove a specific theory, it 1s 
much more difficult to find a theory to fit observed phenomena. 

To return to the initial motivation for basin-scale acoustic tomography, that of ob- 
taining synoptic measurements of basin-scale phenomena, we find ourselves both covering 
new ground and old ideas. Tomographic measurements provide a unique window into 
the understanding of basin-scale processes. Many such processes have been theorized, 
but have no direct observational support. Yet while trying to understand these previ- 
ously unobserved phenomena, we can also obtain a fuller understanding of “well known” 
phenomena (such as tides). 

In this study, we have limited ourselves to a class of interesting observations in a very 
narrow frequency range. The data set however, allows for observations of processes having 
periods from about 1000 hours to 6 minutes. Use of this data in conjunction with point 


measurements in the Pacific is a worthwhile and rich topic of research. 





Appendix A 


Island Trapped Internal Waves 


The solution for the perturbed pressure field of internal waves scattering from a cylinder 
IS, 


Eg au 


Tn o-—oo 


Ja (kr) — ads (ka) — me) Y y (kryl emo, (AD) 
t kal! (ka) — Haka) 7 l 


where the primed variables denote derivatives, Jm and Hm are Bessel and Hankel functions 
of the first kind of order m respectively, A is the amplitude of the incident plane internal 
wave with frequency o and horizontal wavenumber k, a is the cylinder radius, and f is 
the Coriolis frequency (Proudman, 1914; Wunsch, 1972). The origin of the cylindrical 
coordinate system, with radius r and polar angle 0, 1s the center of the cylinder. The 
horizontal wavenumber is related to the frequency through the dispersion relationship for 


internal waves by 
a? = f? 


k* = 
gh 


where h, is the equivalent depth of the n'^ internal wave mode, and q is gravitational 


acceleration. Pseudo-resonances occur if proper choices of k and o result in | P(r,0)| /A 
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greater than about 5 for a given cylinder radius. 

A logical choice for selecting the cylinder radius is made by matching the circumference 
of Oahu to that of a cylinder (29.4 km) (Luther, 1985). Geographical considerations limit 
the maximum island radius to 40 km, at which point the cylinder would touch the 
neighboring island of Molokai. Although not likely, bathymetry suggests the possibility 
of the Oahu/Molokai/Maui group of islands being considered as a single model island, 
resulting in a cylinder radius of about 110 km. 

Historical temperature and salinity profiles near Oahu are used to determine the equiv- 
alent depths of the first three internal wave modes (.88 m, .25 m, and .10 m respectively ). 
We evaluate | P(r,0)| for equivalent depths of .01—1.0 m in .01 m increments and island 
radii of 15-50 km in 5 km increments as well as 90-120 km in 10 km increments. For each 
possible pairing of these equivalent depths and island radii, we evaluate the magnitude 
of the scattering function between 10 and 25 hours at .28 hr increments. The value of 
| P(r,0)| /A does not exceed 3.0 and thus, no significant pseudo-resonances occur. 

We also search for poles of P(r,0) in the complex frequency plane near the real fre- 
quency axis by evaluating a contour integral (Figure A-1) for the same island radii and 
equivalent depth combinations used in the previous calculation. These contour integrals 
are all zero, indicating P(r,0) contains no poles which cause pseudo-resonances at periods 


between 10 and 25 hours. 





Re{a} 
(cph) 





Figure A-1: The contour integral path of P(r,0) encloses a rectangular area in the complex 
frequency domain, where the real axis represents the domain of pure oscillatory waves. 
In order for a pole of P(r,0) to cause a pseudo-resonance at periods between 10 and 25 
hours, it must lie within the region bounded by Q > 5, where Q is the resonance factor 
for a lightly damped oscillator (Q = 5Re{o}/Im{o}). Points on the contour where Q = 5 
are marked with an x. 
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Appendix B 


Reliability of Cross-Spectral 


Phase 


We treat auto- and cross-spectral peaks as isolated events which are too large to be caused 
by measurement error (Spiesberger et al, 1989b). In other words, we are not interested 
in standard calculations of 95% confidence intervals because we do not care if a certain 
acoustic event would be observed in 95 out of 100 different experiments. 

For two sinusoidal signals of the same frequency and specified phase lag (in the absence 
of noise), the cross-spectral phase at the signal frequency is equal to the phase lag between 
the signals. When dealing with real data, however, noise is superimposed causing reduced 
coherency and cross-spectral phase errors. Additionally, some form of averaging must 
be performed on the cross-spectra in order to obtain a coherency spectrum (Jenkins 
and Watts, 1968). This significantly increases the bandwidth of the cross-spectra, with 


the result that cross-spectral phases are biased by noise and other signals within the 





bandwidth of the cross-spectra. 

In light of the above, we estimate the reliability of cross-spectral phases using the 
barotropic tidal signals in the tomographic data. Headrick (1990) has analyzed the 
barotropic tidal signals in the tomographic data, and finds close agreement between the 
data and models of the diurnal tides. A continuous 110-day data segment is available 
from four receivers. Dividing these data sets into consecutive 55-day sections, we com- 
pute cross-spectra for each receiver combination using the same processing and averaging 
used on the 42-day data segments discussed in the body of this paper. At the diurnal tide, 
Ki, cross-spectral phases obtained from each 55-day segment at each receiver combination 
is consistent within +15°. We therefore use 415? as a measure of the reliability of our 


cross-spectral phases. 





Appendix C 


A Two-Mode Gravity Wave 


Model 


Let the observed acoustic travel-time at receiver à be, 


c;sin(ot + 6;), (CHI) 


where c; is the magnitude of the rms travel-time estimated from auto-spectra, 9; is the 
cross-spectral phase (referenced to R3), o is the mode frequency, and t is time. We 
model the acoustic travel-time at receiver 2 from a linear combination of the clockwise 


and counterclockwise travelling resonant barotropic gravity waves (Section 4.4.1) as, 


~ 


¢; sin(o + 6;) = cyan sin(o + $i; + C1) + a2aj2 sin(—o + Qí2 + (2) CESA 


where the variables with tildes are the model estimates, and the riglit side of the equation 
corresponds to the summation of the clockwise and counterclockwise modes. 
Thus for each mode we have four model parameters, (a;,a2, 5 .(2). which are varied 


to minimize the error between the observed travel-times and the model estimates. We 
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establish a relevant (but somewhat arbitrary) cost function, J, which evaluates how well 
the two-mode model fits the observed data at all receivers for a given mode. After some 
trial and error, we selected a cost function which penalizes percentage errors in magnitude 


and deviations in phase between the model and the data. The cost function selected is, 


7 INA 
C; C; 2 2 Yo 
Ha1,02,6,6) = 2 (2 +5) + 100 : (6; — 0;) , (60) 
i=1 1 "t 
where the summation is taken over all receivers. 

For each mode we find, after minimizing the cost function, that the magnitude and 


phase obtained from the two-mode model of resonant barotropic gravity waves performed 


no better than the results using a single mode model. 
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